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испаряющейся жидкости и являющейся моделью высыхания 
невязкой, нетемпературопроводной протяженной капли. Для 
построения автомодельных решений в работе выполнен переход 
к инвариантам Римана. Автомодельные решения представляют 
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собой функции времени и координаты, определяющие высоту регогтеа. Зе-зипи]ае4 зоТанНоп$ аге апсНопз о Яте апа роз!{оп 
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функции, описывающей высоту капли в процессе 
испарения-конденсации. Выявлена область применимости 
различных автомодельных решений к моделированию 
различных ситуаций высыхания капель и пленок. 
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Введение. Математическое моделирование процессов тепло- и массопереноса при испарении-конденсации капель 
жидкостей применимо в настоящее время во множестве современных технологий: медицинской диагностике [1], 
фармакологических исследованиях [2], кристаллографии белка [3], для растягивания ДНК и РНК [4], полиграфии [5], 
создании структурированных поверхностей [6], производстве наноструктур [7], микроматриц, в том числе 
однокристальных лабораторий (1абз-оп-а-сЫр) [8]. В связи с этим в последние годы разработано большое количество 
моделей, описывающих высыхание капель, либо пленок жидкостей. Для математического моделирования процессов в 
испаряющейся капле применяются различные подходы: в первую очередь, те или иные следствия системы уравнений 
Навье-Стокса, а также метод осреднения трехмерных задач по высоте пленки, либо капли испаряющейся жидкости. 

Метод осреднений (аналог уравнений мелкой воды) для моделирования высыхающих капель и пленок использован 
в работах [9-12]. Несмотря на существенные успехи при моделировании конкретных ситуаций в физике, медицине и 
биологии, используемый подход имеет ряд дефектов, таких как необходимость задания плотности пара на поверхности 
капли (на основе эмпирических предположений), потребность во множестве локальных допущений, понижающих уровень 
общности модели, а также преимущественно численный подход к исследованию получающихся начально-краевых задач. 
Процесс осреднения, используемый в перечисленных работах также проводится на основе упрощающих предположений и 
при отсутствии четких математических обоснований. 

В работе [13] разработан универсальный математический аппарат для исследования процессов тепло- и 
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массопереноса в испаряющихся каплях и жидкостях. Введен в рассмотрение оператор осреднения, на основании свойств 
которого осуществлено осреднение и замыкание приближения Обербека-Буссинеска, наиболее часто используемого для 
описания гравитационной конвекции в жидкости. Построен ряд математически корректных моделей высокого уровня 
общности, среди которых наиболее существенной является задача тепловой конвекции, осредненная по тонкому слою 
испаряющейся Жидкости: 


п, +41 (125) = -УФ; (1) 
5, Аза + мА = 0; (2) 
3 3 
Ф, ал + ох = 0; (3) 
с, + 15Ус = Ас. (4) 
Здесь: Й — переменная во времени высота капли, зависящая от координат; 5, ф и с — осредненные по 


толщине капли скорость массопереноса, тепловой поток и концентрация твердой примеси в капле. Величина У — 


безразмерный параметр. Система (1)-(4), описывающая процесс испарения в общем случае, должна быть дополнена 
начально-краевыми условиями, определяемыми из физических допущений для каждого конкретного случая. Начальные 
условия могут определять начальную форму либо высоту капли, а также осредненные значения скорости массопереноса, 
теплового потока и концентрации примеси в начальный момент времени. Выбор краевых условий существенно зависит от 
физических предположений, положенных в основу конкретной модели. В зависимости от степени гладкости основания, 
может иметь место пиннинг границы трехфазного контакта (ее фиксация на шероховатом основании), либо депиннинг (ее 
отрыв на гладком основании — например, слюде, тефлоне). Угол смачивания поверхности будет острым для смачиваемой 
(гидрофильной) поверхности и тупым — для несмачиваемой (гидрофобной); кроме того, во втором случае изначально 
будет иметь место «опрокидывание» профиля капли, что существенно усложняет решение задачи. В случае, когда капля 
примыкает к жесткой стенке, на ней также необходимо задать угол смачивания или (и) условия теплоизолированности и 
непроницаемости для потока жидкости. 

Как отмечается в работе [13], одним из вариантов упрощения процедуры получения решений является 
рассмотрение пространственно-одномерных моделей — для вращательно-симметричной капли или «протяженной 
капли», бесконечно вытянутой в одном из горизонтальных направлений. Оба случая соответствуют общепринятым 
упрощениям при моделировании капель и тонких пленок. Вращательно-симметричная модель капли отвечает случаю 
небольших по размеру капель (то есть радиус капли соразмерим с ее тощиной), принимающих под влиянием сил 
поверхностного натяжения осесимметричную форму. В этом случае в уравнениях модели отсутствует зависимость 
функций от полярного угла и размерность задачи понижается. Это позволяет построить аналитические решения задачи для 
осевого сечения капли. Модель «протяженной» капли соответствует случаю тонкой пленки, обладающей относительно 
большой протяженностью и однородностью в одном из направлений (например, в направлении оси Оу). Предполагается, 


что все физические и геометрические характеристики капли зависят только от переменной х, а решения рассматриваются 
в сечении капли, перпендикулярном оси Оу г Следует отметить, что все краевые эффекты на линии трехфазного контакта 


(жидкость-воздух-твердое основание), возникающие при испарении капли, должны иметь качественное соответствие для 
обоих моделей. 

В работе осуществлено решение задачи (1)-(4) для вращательно-симметричной капли методом характеристик в 
сочетании с численными методами. Установлен ряд интересных эффектов, в том числе зависимость от времени формы 
профиля поверхности капли в окрестности точки трехфазного контакта. Установлено, что в начальный момент времени 
форма профиля близка к линейной; однако с течением времени наблюдается ее искривление, а в некоторый 
фиксированный момент времени происходит «опрокидывание» профиля капли. В работе [13] показано также, что 
сингулярности решений, возникающих при й-»0 для осредненных решений, являются мнимыми и устраняются 
предельным переходом. 

Соответственно, актуальной проблемой является получение аналитических решений для задачи (1)-(4), в первую 
очередь, ее автомодельных решений [14—15]. В соответствии с [14], явление рассматривается как автомодельное, если 
распределения его характеристик в различные моменты времени получаются друг из друга преобразованием подобия. 
Обнаружение автомодельности упрощает вычисления и представление характеристик явления, развивающегося во 
времени. 
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Целью настоящей работы является построение автомодельных решений задачи тепловой конвекции и их 
исследования с точки зрения применимости к моделированию различных ситуаций испарения-конденсации жидкости. Для 
построения автомодельных решений в работе выполнен переход к инвариантам Римана. Получены определяющие 
соотношения, позволившие построить три типа автомодельных замен, понижающих размерность задачи и сводящих ее к 
системе обыкновенных дифференциальных уравнений (ОДУ). Для каждой из систем ОДУ найдены различные 


аналитические решения, определяющие высоту капли й › осреденные по толщине капли скорость массопереноса 5 и 
тепловой поток ф<ф как функции времени г и координаты х. Классификация найденных автомодельных решений 


осуществлена на основании поведения функции Й ‚ описывающей поведение поверхности капли в процессе 
испарения-конденсации. Выявлена область применимость различных автомодельных решений к моделированию 
различных ситуаций, таких как пиннинг (закрепленние границы трехфазного контакта, вызванное шероховатостью 
основания) либо депиннинг (отрыв границы на идеально гладком основании, на котором располагается капля), различная 
геометрическая форма поверхности капли (линейная, выпуклая, вогнутая) и т.п. 

Задача для протяженной волны в инвариантах Римана. Для протяженной капли, то есть бесконечно вытянутой 
в одном из горизонтальных направлений, задача тепловой конвекции (1)-(4), принимает одномерный вид: 


В+ 2И5й, +125, = -№,Ф; (5) 
5, ее =0; (6) 

3155 

2 4 
——Иф5, +—Й5ф, = 0. 7 
а т. о (7) 

Осуществим переход к инвариантам Римана К, Р, О, связанным со старыми переменными соотношениями: 
й= ЕР ЗМ; 5=Р;  ф=ОР. (8) 
В результате система (5)—(7) принимает следующий вид: 

В, +2КР\К, =-ИРО; (9) 
Р + ВР, =0; (10) 

а 
о ве О. =0. (11) 


Для построения автомодельных замен выполним переход к функциям К В Р ь О ь 








в=ок, Р=ЕВР, 0=10, 1=6, — ХЕБХ, (12) 
где а, В, у, 0, = — параметры. При подстановке соотношений (12) в систему (9)-—(11) она приводится кследующему 
виду: 

—_ 2014 м — а 
5 “В Бр“... (13) 
& 
— а 
р 4) 
б 3 
-- дав 
То рб (15) 
б 3 & 
Инвариантность системы (13)—(15) по отношению к (9)-(11) обеспечивается выполнением равенств: 
а _ 2814 а» В _ ав“ . у_ ов!“ 6) 
б 8 | б = б Е” 
На основании (16) получим два определяющих соотношения: 
& 
=; о = Вуб. (17) 
ов! 


Соотношения (17) позволяют получить автомодельные замены, понижающие размерность исходной задачи 
(5)-—(7) и приводящие ее к системам ОДУ. Рассмотрим основные три типа автомодельных замен и соответствующих им 
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автомодельных решений задачи тепловой конвекции. 
Автомодельная замена 1. Полагая © =В =1, из уравнений (17) получим: 


б=$5, =—. (18) 
ие 
Соотношения (18) указывают на существование автомодельных замен задачи: 
Р=,  би=и, (19) 
ГА 
где и и г — новые переменные. Действительно, применение (19) к (9)—(11) сводит систему дифференциальных 
уравнений в частных производных к следующей системе ОДУ: 
К. кр“ = = —ИРи; (20) 
Р| 2 КРМ: |=0: 21 
. з < 1 =У, ( ) 
4 
ини, [+4 ВР 0 (22) 


Система (20)-(22) имеет два решения, определяемых уравнением (21). 
Решение 1.1. Пусть 


3 








ВР =. (23) 
2 
Подстановка (23) в уравнение (22) приводит к решению вида: 
на, (24) 
С; — константа интегрирования, а подстановка (24) в (20) позволяет получить решение системы в инвариантах Римана: 
15 
81 5 
К=| С. -— СМ 3 25 
2—3 о (25) 
—45 
81 4 81 5 
Р=— < | © -—С М } 26 
16 230 о (26) 
5 
о=с-. (27) 
Вернемся от инвариантов к обычным переменным, и на основании (25)-—(27) получим решение исходной системы (5)-(7): 
8 81 4/5 
=— [с си" . (28) 
2х 16 
—45 
81 4 5 81 5 
= СЕ -— СУ х . 29 
16 Е (29) 
81 х 5 81 5 в 
=—С СУ СУух : 30 
ом (©, те "т" (30) 


Функция й ‚› определяющая тощину капли, ведет себя по-разному в зависимости от знаков констант С ы С. . 


Пусть С.С. >0 (\>0). Тогда из (28) видно, что функция Йй обращается в ноль при 


ее ЗО (31) 
3\|2,С, 


Точка (31) соответствует краю капли, причем это значение увеличивается с возрастанием Е. Из формулы (28) 





следует, что при С.С. >0, если О<х<х, ‚ то й, <0, й/>0 то есть имеет место убывание функции по х и 
возрастание г (конденсация). Если х>х,, то й’ >0, Й, <0, что означает возрастание функции по х и убывание 1 


(испарение). Таким образом, решение (28)—(30) может быть использовано непосредственно для описания испаряющейся 
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капли, если в качестве начального момента времени положить {= 1 > 0. Предполагается, что левый край капли (стенка, к 
которой примыкает капля) находится в точке х = ху > х,. Тогда константы С:, С› можно определить, например, задавая 
значения И(хо) =, 5(х0)=50 при 1=Ю. Если в качестве левого края капли взять точку х=0, то формулы (28)-(30) 


определяют конденсацию капли. Однако их нельзя использовать целиком из-за сингулярности. Используем сшивку 
полученного решения с константой (рис.1): 


Вх) 





х К х К хо х 


Рис.1. Поверхность капли, описываемая Решением 1.1, С.С >0 (У >0). 
хх — точка трехфазного контакта. Левая область — конденсация капли; х =0 (начало кооординат) соответствует жесткой 
стенке, к которой примыкает капля. Правая область — испарение; х = ху — жесткая стенка. 


Кривые (1) и (2) соответствуют различным моментам времени 1 <; 





4/5 
1 
. [се-1 см» 5 О<х<х,, 








- _ 
в. 27% 16 (32) 
8 $1 . 4/5 
273 ср вт г Хс <х< Хх. 
—4/5 
вх су" си Фа. 
_ (33) 
О: г 
2" С» в ; Хх. <х< ху. 
] 5 ] —4/5 
а с бо 
16 1Р Ю^” : 
а (34) 


х 81 Ра 
ен. , ое. 
= 2 16° с [3 


Таким образом, при О<х<х, функции й,5,ф представляют собой константы, зависящие от времени. 
Соотношения между константами, входящими в (32)-(34) можно получить, задавая й(0,1) =, 5(0,0)=50 , 
ф(0,5) = фо. Следовательно, Решение [1.1 при С!,С. >0, определяемое формулами (32)-(34), в различных областях 


задания может быть использовано для описания как конденсирующейся, таки испаряющейся капли. 
Если же С, С. <0 › ТО, области при х<х, и х>х, меняются ролями: левая часть отвечает за испарение, а 


правая — за конденсацию. Если С, >0, С, <0 (либо С, >0, С, <0), то подкоренная функция в ноль не обращается. 


2_|9С> 
Однако у функции Й имеется минимум в точке хи =-—5|-— = 
32С, 


случае описывает испарение (конденсацию) пленки жидкости, расположенной неравномерным слоем между двумя 


1. Можно предположить, что функция Й в этом 
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твердыми стенками. 


Решение 1.2. Рассмотрим второе решение уравнения (21), соответствующее Р. =0: 


Р=р,, (35) 

где О, — константа интегрирования. Тогда оставшиеся уравнения системы преобразуются к виду: 
В. (2вр, -=)=-И Буи; (36) 
ини. [2+8 |=0 (37) 


Система (36)-(37) допускает построение асимптотического решения в виде рядов. Возвращение к исходным переменным 
приводит к квадратичной асимптотике автомодельного решение системы (5)-—(7): 








2 
20 р, В 
й—| В+ |= дов +0) : В, =58 +2М.Бищ. (38) 
1 1 
$=ра, (39) 
4 4 2 
р 
а ТЫ (40) 
1 460 Е 8 В 


Исходя из требования неотрицательности функции толщины капли Й в начальный момент времени, будем 
считать, что ОА, >0. Предположим для определенности: О; >0, К >0. Парабола (38) имеет вершину в точке 
Кок, +2\ БЕ 

ИБ Чо 





в(- ДКК; Ко + 2% к, }, и пересечение с осью ОХ вточке х, =ИК| - БВ + 


Таким образом, с течением времени высота капли, описываемой (38)-(40) возрастает, а площадь ее увеличивается. 
Следовательно, автомодельное Решение 1.2 соответствует конденсации капли (Рис.2). Неизвестные константы ид и Ку 


могут быть найдены, например, из задания углов смачивания в начальный момент времени [= Ю >0 точках х=0 и 


хХ=хь: №, (х=0) =1804, й(х=хь) = ва. 





|) .х_К хк х 
Рис.2. Поверхность капли, описываемая Решением 1.2, конденсация, ПД, >0, КЮ >0; х= 0 — жесткая стенка, 
к которой примыкает капля, х = х, — точка трехфазного контакта. 


Автомодельная замена 2. Пусть В=у=1. Тогда из (17) следует, что 


5=а —5=5°. (41) 
Следовательно, на основании (41) применимы автомодельные замены: 
3% 
==,  К=И, (42) 


1 
приводящие систему (9)—(11) к виду: 


У+2У, (-=+УРМ) = М РО; (43) 
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2 - 2: дур") = (44) 
о.- 22+ чур") =0. (45) 
Из уравнения (44) следует, что 
Р. =0, (46) 
ИЛИ 
УР, (47) 
а из уравнения (45) — что 
0. =0, (48) 
либо 
ур” 52 (49) 


Рассмотрим решения, соответствующие различным комбинациям (46) - (49). 
Решение 2.1. Пусть выполняются (46) и (48), следовательно, 


Р=Р', (50) 
О=О,, (51) 
гле Р, О, — константы. Тогда решение системы (43)-(45) имеет вид: 
В 1 . 
х+?| УР?О, +1 | = РУ| ВР +1 У,В*0,+— ||, (52) 
и У 
$=Р,, (53) 
0=0:В", (54) 
где Р, О, У — константы. Кривая (52) есть парабола с вершиной в точке 
РО, +1 
В -пвк; № | кет (55) 
В у 


2рб .. 
и пересечением с осью ОХ в точке С(Е”Р’КУ,О;;0). Если, например, РК <0, то вершина расположена в первой 
координатной четверти, а ветви параболы повернуты влево. Таким образом, в соответствии с (55), в начальный момент 
времени {=0 получаем параболу с вершиной в начале координат: х=Р'УЙ?. При увеличении времени вершина 


парболы смещается вправо и вверх. Таким образом, наблюдается опрокидывание профиля капли, а затем высота капли 
увеличивается. Следовательно, Решение 2.1, определяемое формулами (52)-(54), соответствует конденсации капли на 
вертикальной стенке, совпадающей с осью ординат (рис. 3). 


|9) 





х 


Рис. 3. Поверхность капли, описываемая Решением 2.1, конденсация капли на вертикальной стенке. 


Линии соответствуют различным моментам времени Ги] < ь 
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3 
Решение 2.2. Пусть теперь выполняются (46) и (49). Следовательно, Р = РА, У = Ш где РА — константа. В 


1 
этом случае получим второе решение системы (43)—(45): 
Пи — 28%, ба — 

221 УВЕ 


Решение 2.2. определяемое уравнениями (56) соответствует случаю испарения, когда поверхность капли линейна, 








(56) 


жидкость уменьшается со временем. В данном решении по смыслу задачи РО, >0; в противном случае формула (56) 


описывает «перевернутую» каплю. 











4 
Решение 2.3. Пусть теперь выполняются (47) и (48), О=0О; Р= ыы . В этом случае авнение (43) 
У р 1 уз учае ур 
приобретает линейную форму: 
ау У 810,53 
+—=-И : 57 
@ 4  ° 4% ее 
Решив (57), получим третье решение системы (43)-—(45): 
14 27 эм * 
В=к Ч, 4 = 5 А: И , (58) 
Р=81574 ", (59) 
О=О0.. (60) 
После возвращения к исходным переменным (8) в (58)—(59), имеем автомодельное решение исходной задачи: 
15 хИ5 57 4/5 
де = 5 [и -Эбикм | (61) 
1х4 27 —45 
Е 8 а 5-45 д о , (62) 
1 7 
=. 27 — 
Ф=О, 5-45 А? Ау ом | (63) 
#45 7 


Функция й, определяемая формулой (61) не имеет сингулярностей ни по х,ни по г. При фиксированном 1 


421 
ТА! 2 
функция имеет пересечения с осью ОХ вточках х„=0, и х, = 7.0, 1“ и максимум в точке 
ое 
ВА 421 
а , (О<х, <, ). 
675% О! 


в(х) 





Рис.4. Поверхность капли, описываемая Решением 2.3. х, — точка трехфазного контакта. 
Левая область — конденсация капли; х = 0 — фиксированный левый край капли; х = х; — движущийся правый край. 
Правая область — испарение; х = х) — жесткая стенка, к которой примыкает капля. Кривые (1) и (2) 


соответствуют различным моментам времени Ё < 15 
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Таким образом, решение (61)-(63) может быть использовано в области х<х<олх, для описания 
конденсирующейся капли в условиях пиннинга ее левого края, а в области х, <х<ху — испаряющейся капли (ху — 


жесткая стенка, к которой примыкает капля) (рис.7). Соотношения между константами, входящими в (61)-(63) можно 
получить, задавая й(х„, 0) =, 5(хт, 0) = 50. Следовательно, Решение 2.3 в различных областях задания может быть 


использовано для описания как конденсирующейся, таки испаряющейся капли. 


Автомодельная замена 3. Пусть © =у=1. Тогда из (17) следует, что 


ВТ, веб, (64) 
б 
Следовательно, на основании (64) применимы автомодельные замены: 
х 
<= Зи’ Р=\/, (65) 
где с и \ — новые неизвестные переменные. Применение формул (65) к системе (13)-(15) приводит ее виду: 
®.[- = =+ али = ИМО; (66) 
ин. [Че+акии 0 (67) 
4 3 
3 4 рим 
——&+—А\ |=0. 68 
Из (68) следует, что 
О. =0, (69) 
или 
щ_ 9 
В" =— 2. (70) 
16 


Система (66)—(68) имеет два решения, определяемых уравнениями (69) и (70). 
Решение 3.1. Рассмотрим вначале случай (70). Система сводится к виду: 








в ЗЕМ; (71) 
3 
и’ [3 =. (12) 
8 
Решение уравнения (72) имеет вид: 
И =В, ечо(-4 2? (73) 
В, — константа интегрирования. Тогда из уравнений (70), (71) и (73) получим: 
9: 1 2 
К = ехр| —<” |; 74 
16! 23 у 
275 2 › 52 
@= [= с }= (5: } (75) 
маи 3 Е 
Вернемся в (74)-(75) от инвариантов к обычным переменным (8), получим: 
9х 4х? 
й=— еф ——, 76 
В. \ зу _ 





2, 
‚в Е | (7) 
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275 О" Г 
ф= 1+ ехр р (78) 
А 3/3 |=; 


Из формул (76)—(78) следует, что функция й возрастает по х, убывает по #; функция 5 убывает по х, убывает по 1; 





функция ф возрастает по х, убывает по 1 (рис.5). Константа В, может быть найдена из задания, например, угла 


смачивания при х= 0. 


В(х) 
1 








о х 


Рис. 5. Поверхность капли, описываемая Решением 3.1, конденсация капли в условиях пиннинга. 


Линии соответствуют различным моментам времени |8 < 15 


Таким образом, Решение 3.1. соответствует испарению в условиях пиннинга, то есть неподвижной границы трехфазного 
контакта в точке х=0. 
Решение 3.2. Наконец, на основании (69) получим частное решение задачи: 


4/5 4/5 

3 х 63 х 
Ех, = В, ф=-— в, (79) 
В 1 З\В + 
где В есть произвольная постоянная. Решение 3.2, определяемое (79), соответствует конденсации капли в условиях 
пиннинга (рис.5). 


Заключение. Как следует из представленного материала, задача тепловой конвекции обладает большим 
количеством автомодельных решений, получаемых при переходе к инвариантам Римана. С ПОМОЩЬЮ подбора констант, 
соответствующих краевым условиям, данные автомодельные решения можно применять для моделирования различных 
практических ситуаций, таких как испарение капли и ее конденсация в различных геометрических конфигурациях: 
пиннинг-депиннинг (в том числе только на части границы); задание различных форм поверхности капли и начальных точек 
трехфазного контакта, задания начального угла смачивания и т.д. 
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